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ABSTRACT 



Aims. In the first few minutes of a newborn neutron star's life the Hall effect and differential rotation may both be important. We 
demonstrate that these two ingredients are sufficient for generating a 'shear-Hall instability' and for studying its excitation conditions, 
growth rates, and characteristic magnetic field patterns. 

Methods. We numerically solve the induction equation in a spherical shell, with a kinematically prescribed differential rotation 
profile Cl(s), where s is the cylindrical radius. The Hall term is linearized about an imposed uniform axial field. The linear stability of 
individual azimuthal modes, both axisymmetric and non-axisymmetric, is then investigated. 

Results. For the shear-Hall instability to occur, the axial field must be parallel to the rotation axis if Q.(s) decreases outward, whereas 
if Q.(s) increases outward it must be anti-parallel. The instability draws its energy from the differential rotation, and occurs on the 
short rotational timescale rather than on the much longer Hall timescale. It operates most efficiently if the Hall time is comparable 
to the diffusion time. Depending on the precise field strengths Bo, either axisymmetric or non-axisymmetric modes may be the most 
unstable. 

Conclusions. Even if the differential rotation in newborn neutron stars is quenched within minutes, the shear-Hall instability may 
nevertheless amplify any seed magnetic fields by many orders of magnitude. 
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The classical dynamo mechanism for generating cosmic mag- 
, netic fields is based on the inductive effects of rotating 
. anisotropic turbulence. Because the same turbulence prevents 
uniform rotation from being a solution of the Reynolds equation, 
in the majority of dynamo models, the basic mechanism for pro- 
ducing toroidal fields from poloidal fields is differential rotation. 
However, it alone only produces toroidal fields that grow linearly 
. in time, until the resulting Maxwell stresses suppress the differ- 
ential rotation. To produce a self-excitation of magnetic fields as 
an instability that grows exponentially in time, another effect is 
. needed to produce poloidal fields from toroidal ones. 

There are many possibilities. The simplest one is a merid- 
ional flow. For certain combinations of differential rotation, 
meridional flows, and finite conductivity, the coupled system al- 
, lows the self-excitation of (purely non-axisymmetric) fields (e.g. 
iDudlev & Jamesiri989h . It is also possible to close the circuit by 
the a-effect, which represents the existence of an electromotive 
force parallel to the mean magnetic field. The result is a so-called 
af2-dynamo, which for a high magnetic Pra ndtl number can 
be modified by the actio n of meridional flow (IChoudhuri et al.l 
[T995HKiiker et al.ll2001h . Driven turbulence subject to a shear 
flow in liquid metals is also able to work as a dynamo if the tur- 
bulence intensity is stratified in the direction o rthogonal to the 
shear-flow plane dRiidiger & Kitchatinovll2005h . 

The dynamo equations resulting from purely hydrodynamic 
flows are all linear, and the amplification can start from arbitrar- 
ily weak seed fields. This is not the case, however, for possi- 
ble dynamo mechanisms where the coupling Btor — > Spoi de- 



pends on magnetic instabilities. Such alternative mechanisms 
only work if the seed field is sufficiently strong. There are a va- 
riety of magne tic instabilitie s on wh ich the dynamo mechanism 
can be based. ISchmitt et al.l (1199 6*) started the construction of 
such models with their on-off dynamo model where the a-effect 
only exists if the magnetic field exceeds a certain threshold value 
determined by the magnetic buoyancy instability. Other possibil- 
ities have been suggested on the basis of both the magnetorota - 
tional instability (MRI) and the Tayler instability ('S pruitll2002l) . 

In this Letter it is demonstrated that in the presence of a shear 
flow, including the Hall term in Ohm's law is enough to pro- 
duce a 'shear-Hall instability' (SHI). The Hall effect produces 
off-diagonal elements in the conductivity tensor that are of the 
first order in the magnetic field. The SHI is thus a magnetically 
nonlinear instability, but with the particularity that there is no 
threshold value. That is, if the shear is strong enough, the Hall 
effect, hence the initial seed field, can be arbitrarily weak. 

The natural application of this shear-Hall instability is in 
newly born neutron stars. In the first few minutes of its life, 
the crust of the star can rotate differentially, as shown in 
a number of numerical comp u tations of the supernova core 
yMonchmever & Muller '1989'; 'Janka & Moenchmever' '19891 
Dimmelmeier et al. 2002; Kotake et al. 2004; Ardeljan et aTj 
120051 ; IBurrows et alj|2007l) . Differential rotation in the core of 
the collapsing progenitor will ind uce strong diffe rential rotation 
in the newly formed neutron star (lOtt et al.l2005h . 

Initially, the star cools down by neutrino emission. As a re- 
sult of this cooling, the core becomes distinct from the outer 
layer, and the crust co ndenses at a temperature of a few 10"^ K 
dDall'Osso et al.ll2009l and references therein). Any differential 
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rotation can thus only exist in a limited time-interval that de- 
pends on the type of reactions causing the cooling. In the case of 
the direct URCA process, the cooling is extremely fa st, with es- 
timate s giving a characteristic time of a few minutes dPage et alJ 
l2006h . 

However, the very rapid rotation rates, with periods of order 
10"^ s, mean that even just a few minutes is long enough for the 
star to undergo several thousand revolutions. Furthermore, we 
show that the SHI grows on this very fast rotational timescale, 
so even though it only exists for a few minutes, this is sufficient 
for it to amplify a seed field by many orders of magnitude. Of 
course, the time evolution of the shear is excluded in this simpli- 
fied model, so the exact amplitude of the field at the end of the 
SHI phase must still be deduced from full nonlinear simulations 
of the Hall-MHD system. 



2. The shear-Hall instability 

The Hall effect is important whenever a plasma is magnetized 
enough. The level of magnetization is determined by the Hall 
parameter, Rb - crBo/nec, with Bq a characteristic magnitude 
of the field, cr the electric conductivity, n the electron number 
density, e the electron charge, and c the speed of light. If the 
Hall parameter Rb is large, the Hall effect is important. 

In a neutron star, both the electron number density and the 
conductivity vary with depth. The conductivity also varies with 
temperature. To focus on the basic physics of the instability pro- 
cess, we neglect all such variations and treat both n and cr as con- 
stants. However variations in n in particular can i nteract with the 
Hall effect in interesting and unexpected ways ( iVainshtein et alj 
120001: iHollerbach & Riidig"en2004.) . 

The magnetic Reynolds number is Rm - ^o/t], where R 
is the star's radius, Qq the maximal angular velocity, and rj its 
magnetic diffusivity. Scaling length by R, time by R^/rj, u by 
RQ.(i, and B by Bq, the induction equation becomes 

dt ~ 

-RbV X ((V X B) X B) - V X (V X B) -h RmV x (m x B). (1) 

We linearize this equation about the uniform axial field Bq - 
e,, yielding 

dt ~ 

-RbV X ((V X B) X e^) - V X (V X B) -H RmV x (« x B). (2) 



The flow u is kinematicaUy prescribed as 



u - sQeA 



where Q 



(3) 



where s is cylindrical radius, and sq = 0.5. By itself such a flow 
generates no interaction with the imposed field Bq - e,, not 
even a linearly growing toroidal field. It is only with the addition 
of the Hall term that instabilities may develop in this otherwise 
stationary basic state. 

Investigations of the full fi eld/flo w cou pling have already 
been conducted by IWardld d 19991) and iBalbus & TerquemI 
d200lh in the loc al approximation, as well as globally by 
iRiidiger & Kitchat inov (2005). These works show that the inter- 
action between the Hall effect and the MRJ may lead to both sta- 
bilization or destabilization of the system, depending on the val- 
ues of the Lundquist number (S - BqR/ -\/(47Tp)ti, where p is the 
mass density) and Rb. In the context of IRiidiger & Kitchatinov! 



d2005h . our case corresponds to the limit of small S. For exam- 
ple. Fig. [T]illustrates how the increase in S (the influence of the 
Lorentz force) has a stabilizing effect. 

It is also useful to derive the energy equation associated with 
either Eq. ([T) or (|2]i. Taking the dot product of either equation 
with B and integrating over the volume, one obtains 



J ^B^dV = - Jc^xBfdV +Rm J B-Vx{uxB)dV. (4) 

The Hall term exactly conserves magnetic energy, in both the 
fully nonlinear and the linearized versions of the induction equa- 
tion. The only source of energy is from the differential rotation. 
The Hall term is therefore essentially a catalyst, necessary for 
the instability to proceed, but energetically incapable of driving 
instabilities itself ( Warein g & HoUerb ach 2009). This situation 
is exactly analogous to the MRJ, where the energy source is also 
purely from the differential rotation, but a magnetic field is nev- 
ertheless necessary. It is, perhaps, not surprising then that the 
shear-Hall instability, like the MRI, grows on the fast rotational 
timescale, as we demonstrate below. 

According to the local analysis of lUrpin & Riidige^ d2005h . 
for a general differential rotation depending on both s and z, a 
necessary condition for the existence of unstable modes is 

{k-B)\k,— -k,— \<Q, (5) 

where k is the wavenumber For Q - Q.{s), this reduces to 

B.^<0, (6) 
as 

so that, for a given gradient of the angular velocity, instability 
only occurs for one sign of (Wardle 1999; Balbus & Terquem| 
2001). For our particular profile (O this condition is only ful- 
filled for positive B~, hence Rb. For negative Rb there is no SHI 
instability. 

Because Eq. (|2]i is linear with an axisymmetric basic state, 
the different azimuthal modes e\p(im<p) decouple and can be 
studied separately. Each azimuthal mode further decouples into 
equatorially symmetric and anti-symmetric modes. We will 
present results for SO, AO, SI, and Al, where 'S' or 'A' refers 
to the radial component of the field being equatorially symmet- 
ric or anti-symmetric, and the '0' or ' 1 ' refers to the wavenumber 
m. 

Equation d2|l is solved in a spherical shell using the spher- 
ical harmonics code described by HoUerbach (2000). At both 
Rin = 0.7 and /?out = 1 boundaries, the solution smoothly 
matches a potential external field. At /?out this is indeed the cor- 
rect condition, but at it is clearly a considerable simplifica- 
tion of the true physics governing the coupling to an internal 
core field. Resolutions up to 30 Chebyshev polynomials in r and 
70 Legendre functions in were used, which was sufficient to 
resolve the eigenmodes up to Rb = 8. The Hall parameter close 
to unity is appropriate during the initial stages of a neutron star, 
and also turns out to be where the shear-Hall instability operates 
most efficiently. 

3. Results 

3.1. Stability maps 

Figure [1] shows the instability curves for the four modes SO, 
AO, SI, and Al. Over most of the Rb range the axisymmetric 



T. Kondic, G. Riidiger and R. HoUerbach: The shear-Hall instability 



3 



modes are slightly preferred over the non-axisymmetric ones, 
but the curves are otherwise all qualitatively similar. The min- 
imum value of Rm occurs when Rb - 2 - 4. At this minimum 
the shear frequency lOsh is arou nd 400 times the Hall fr equency 
a)}i, consistent with the result of lUrpin & Riidige"3(l2005l) that the 
instability only exists for l^shl > I'^'hI- 

If Rb is smaller than ^ 2 - 4 the system tends toward the (sta- 
ble) pure differential rotation case, thus explaining the asymp- 
tote Rm ^ oo as Rb — > 0. At the other end, as Rb becomes 
large the Hall effect increasingly dominates, but as noted above, 
it is energetically incapable of driving any instability itself, so 
Rm must again increase to compensate. For both Rb 1 and 
Rb » 1, therefore, the differential rotation must become increas- 
ingly strong for the instability to operate. 
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Fig. 1. Stability curves for the axisymmetric (SO and AO) and 
non-axisymmetric (SI and Al) modes. The single points belong 
to the system where the flow equation was treated in addition to 
the induction equation. 



An intriguing feature of the small Rb ends of all four insta- 
bility curves is that none of them shows any evidence of tilting 
back over to the right (down to our lowest investigated Rb value 
of 0.4), that is, there is no indication of the instability switching 
off again for sufficiently large Rm. This is perhaps not so sur- 
prising for the axisymmetric modes; For example, the classical, 
axisymmetric MRI has the same feature. For nonaxisymmetric 
modes, though, it is often the case that a s trong enough differ- 
ential rotation suppresses them again (e.g. Riidige r et al.ll2007t 
^ollerbach et al. 2010). It is not entirely clear why this does not 
occur here - or perhaps it does, but only at much higher values 
of Rm than what we were able to compute here. 

It is also interesting that there is no region where a single 
mode dominates: both A and S modes, and both axisymmetric 
and non-axisymmetric ones, are excited at comparable parame- 
ter values. This suggests that the fully nonlinear regime might 
consist of a mixture of several modes, which can also be consid- 
ered as an 'oblique rotator' . If different modes are indeed present 
simultaneously, the field configuration and amplitude in the two 
hemispheres would necessarily be different. If the heat trans- 
port is furthermore dominated b y the magnetic suppression of 
the thermal conductivity tensor ( Schaaf] ! 19901: iHevl & HernquistI 
I2OOI), then the Hall effect would th ereby lead to different te m- 
peratures, and hence X-ray emission (iBecker & Pavlovl2002h . in 
the two hemispheres. 
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Fig. 2. 7 is a nondimensionalized growth rate (7 = Tioty, where 
y is in physical units), of the SO and S 1 modes, for Rb = 1 and 5. 
The profiles for AO and Al are very similar, and are not shown 
here. 



3.2. Growfh rafes 

To assess the possible astrophysical significance of the shear- 
Hall instability, we need to know not only for what parame- 
ter values it exists, but also what its associated growth rates 
are. Figure |2] shows the growth rate scaled on the rotational 
timescale. (That is, if ydiji is the growth rate on the chosen dif- 
fusive timescale R^lr], then y - (2;7r/Rm)ydiif is the growth rate 
on the rotational timescale Tiot = 27r/Qo.) For both Rb = 1 and 
Rb = 5, for large enough Rm these growth rates y eventually 
saturate at 0(1) values. As before, for the initial onset, all four 
modes behave comparably, although the axisymmetric modes do 
have higher growth rates. There is no preferred equatorial sym- 
metry, because the S and A modes are very similar Surprisingly, 
even though the stability map shows that a nonaxisymmetric 
mode S 1 is less stable in the Rb region from 4 onward, its growth 
is only stronger near the stability line. For higher Rm, the ax- 
isymmetric modes are more dominant. 

Of the three timescales that are 'obviously' present in the 
problem, namely the fast rotational timescale Tiot, the much 
slower diffusive Tjiff = /?^/'7 = (Rm/27r)rrot, and Hall Teaii = 
T"diif/RB timescales, the shear-Hall instability thus grows the 
fastest, and is independent of both Tdiff and THaii. This is sur- 
prising, given that it does clearly depend on both the shear-effect 
and the Hall effect. On the other hand, as noted above, its en- 
ergy source is entirely drawn from the differential rotation. In 
this regard it is quite different from a classical auj dynamo, for 
example, where the growth rate is the geometric mean of Trot and 
Ta, but both the differential rotation and a also contribute to the 
energy balance. 



3.3. l/l/ave numbers 

The spatial properties of the eigensolutions are given in Fig. 
[3] The color map represents the toroidal field structure, while 
the contours depict the field lines of the poloidal component. 
Comparing the two plots, at Rm - 2000 and 10000, suggests a 
transition to smaller scales as Rm increases. The wavenumber 
also increases for (turbulent) mean-field dynamos at increasingly 
supercritical dynamo numbers. However, thanks to the high val- 
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ues of their eddy diffusivity, turbulent dynamos operate with 
moderate values of Rm. The SHI though may easily operate with 
very high Rm, and may play a dominant role in driving MHD 
turbulence if Rm » 1 . 

The relative amplitudes of the toroidal field never exceed 
that of the poloidal field by more than one order of magnitude. 
Independent of the ratio of the Hall time and the rotation time, 
the coupling between the components via the Hall term enables 
their simultaneous, and comparable, amplification. 




Fig. 3. The axisymmetric solution AO for Rb = 1, and 
Rm = 2000 on the left and Rm = 10000 on the right. The other 
modes have the same pattern of exhibiting finer structure in z for 
larger Rm. 



4. Discussion 

It remains to be seen to what extent the SHI is important in 
the regime of newborn neutron stars. We first need to establish 
the reasonable values for the Hall parameter and the magnetic 
Reynolds number. It turns out that owing to a low diffusivit y 
iri = IQ-^^ -^f cm^ s"') (as treated in the code of'Potekhin"2008'), 
the microphysical Rm is orders of magnitude above the stabil- 
ity line in Fi g. [T| If the efi'ec ts of turbulence are important, ac- 
cording to Naso e t all (120081) diff^usivity is much larger and the 
corresponding Rm is lower (~ 10^), but still supercritical. 

The onset of the shear-Hall instability requires a seed field 
of the magnitude determined by the actual Rm. This field can 
be provided by various other proposed mechanisms, fro m flux- 
conservation to dynamo processes (for an overview, see ISpruitI 
[2009). The turbulent Rm would require Rb > 1 , which, for the 
temperatures of lO^ - '" K and densities of lo'" ■ g cm"-', corre- 
sponds to B > 10'^ G. If no turbulence is present, the seed field 
can be much smaller provided that the stability line remains uni- 
valued when Rb tends to zero. 

In the situation where the magnetic fields are strong, there is 
a coupling between the fields and the differentially rotating flow, 
which leads to a stable state of uniform rotation. The timescale 



for th is process is ta - RyfAnplB, the Alfven time (IShapirol 
which is only ~ 10 s for the neutron star conditions. Is the 
SHI quick enough to be of any significance? Taking the fast rota- 
tion of neutron stars into account , our results (Fig. |2]i show that 
the field can easily amplify several orders of magnitude within 
the constraining time interval. 

Because we are convinced that the instability conditions are 
fulfilled and the growth is quick enough, there is only the field 
geometry left to analyze in the context of the standard, inclined 
magnetic dipole model, which is used to explain the spin-down 
of pulsars. The SHI generates, as shown, small-scale fields for 
large Rm. The presence of turbulence is, it seems, a necessary re- 
quirement for producing a global field via this instability. On the 
other hand. X-ray spectra of pulsars sometimes indicate an addi- 
tional, strong small-scale surface component (iPavlov et al.l2002t 
Sanwal et al. 2002). It is possible that the SHI is involved in gen- 
erating the small-scale structure, while the large-scale compo- 
nent should be attributed to another process. 

Rapid growth, as well as the easy excitation, indicates that 
the shear-Hall instability may play an important role in the evo- 
\ lution of the field of newborn neutron stars, prompting further, 
I nonlinear analysis of the field/flow coupling and raising ques- 
I tions about its interaction with other field-generating mecha- 
[ nisms. 
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